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ABSTRACT 

The main observed properties of Ly-a absorbers are investigated on the basis of the- 
oretical model of formation and evolution of DM structure elements. This model is 
generally consistent with simulations of absorbers formation and with statistical de- 
scription of structure evolution based on the Zel'dovich theory. The analysis of redshift 
variations of comoving linear number density of absorbers was performed in our pre- 
vious paper. 

We show that the observed characteristics of Doppler parameter can be related to 
the size of DM structure elements what allows us to explain the observed distribution 
of Doppler parameter. This distribution is found to be consistent with the Gaussian 
initial perturbations. The observed characteristics of entropy and column density, Nui, 
confirm that merging of pancakes is the main evolutionary process at redshifts z > 2. 
The observed sample of absorbers characterizes mainly the matter distribution within 
large low density regions and therefore it is difficult to reconstruct the density field 
from the distribution of absorbers. 
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1 INTRODUCTION 



The intergalactic nature of Ly-a forest was established by 
Sargent et al. (1980) and later many models of absorbers for- 
mation and evolution were proposed (see, e.g., Rees 1986, 
1995; Ikeuchi & Ostriker 1986; Bond, Szalay & Silk 1988; 
Miralda-Escude et al. 1996; Hui et al. 1997; Nath 1997; 
Valageas, Schaeffer, & Silk 1999). The essential progress was 
reached through simulations of dynamical and thermal evo- 
lution of gaseous component with the CDM-like power spec- 
trum which is probably responsible for the formation of ob- 
served galaxy distribution. Such simulations (Petitjean et al. 
1995; Hernquist et al. 1996; Bond & Wadsley 1997; Zhang et 
al. 1997, 1998; Theuns et al. 1998, 1999; Bryan et al. 1999; 
Dave et al. 1999; Weinberg et al. 1998; Machacek et al. 2000) 
reproduce successfully the main observed properties of ab- 
sorbers and connect this problem with a more general prob- 
lem, that is the nonlinear evolution of initial perturbations 
and formation of large scale matter and observed galaxy dis- 
tribution. This progress allows us to consider the properties 
of absorbers in the context of nonlinear theory of gravita- 
tional instability and statistical description of formation and 
evolution of DM structure (Zel'dovich 1970; Shandarin & 



Zel'dovich 1988; Demianski & Doroshkevich 1999, hereafter 
DD99; Demianski et al. 2000, hereafter DDMT). 

It is commonly recognized that in the DM dominated 
universe the fundamental properties of observed matter dis- 
tribution are determined by the evolution of DM struc- 
ture elements. The approximate statistical description of 
expected characteristics of DM structure elements based on 
the Zel'dovich theory was given in DD99 and DDMT for 
the CDM-like transfer function (Bardeen et al. 1986) and 
the Harrison - Zel'dovich initial power spectrum. This pro- 
cess can be outlined as a random formation and merging of 
Zel'dovich pancakes, their successful transformation to fila- 
mentary component of structure, and the hierarchical merg- 
ing of both pancakes and filaments to form rich walls. All 
steps of this evolution are driven by the initial power spec- 
trum. 

Some results of this statistical description were com- 
pared with simulations at small redshifts (DD99, Doroshke- 
vich et al., 1999, hereafter DMRT; DDMT). This comparison 
shows that the main simulated and observed structure char- 
acteristics are consistent with the theoretical expectations. 
Here we use this approach for the analysis and interpreta- 
tion of the absorbers observed at large redshifts as a Ly— a 
forest, and we show that some of the observed characteris- 
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tics of Ly-a absorbers can also be successfully described in 
the framework of this theoretical model. 

The validity of this description, when it is applied to 
structure at high redshifts, is not yet reliably verified with 
available simulations due to a small density contrast of poor 
structure elements dominating at higher redshifts. The first 
tests revealed, however, the existence of three kinds of struc- 
ture elements, namely, high density filaments, clumps, and 
low density pancakes in DM spatial distribution at z = 3. 
The theoretical description (DD99) confirms also the self- 
similar character of structure evolution (at least when the 
Zel'dovich approximation can be applied). More detailed 
statistical comparison of absorbers characteristics simulated 
at high redshifts with theoretical expectations is however 
required. 

Such approach implies more traditional investigation 
of properties of discrete absorbers rather then a continu- 
ous non-linear line-of-sight density field (see, e.g., Weinberg 
ct al. 1998) and, so, we will concentrate more attention on 
statistics of discrete absorbers. The observations of galaxies 
at large redshifts verifies that strong nonlinear compression 
of DM and gaseous components occurs even at z « 3 - 5 and 
earlier (Steidel et al. 1998; Dey & Chaffee 1998; Fan et al. 
2000) and, so, discrete high density absorbers can already 
exist at such redshifts. This approach allows us to reach 
more clarity in the description of formation and evolution of 
structure and to establish correlations between properties of 
observed absorbers and invisible DM component. 

It can be expected that absorbers are predominantly 
related to the more numerous population of moderately rich 
pancakes and to periphery of richer pancakes and filaments. 
This means that information obtained from such traditional 
methods of analysis is related with the extended lower den- 
sity regions of the universe rather then with the rich wall-like 
condensations. In this respect, the absorbers characteristics 
are complementary to the information obtained from the 
analysis of large scale galaxy distribution. In order to dis- 
criminate between properties of absorbers that can be as- 
sociated with the evolution of DM component and are spe- 
cific of the gaseous component it is necessary to compare 
the expected characteristics of DM structure and observed 
discrete Ly-a absorbers. Thus, some fraction of weaker ab- 
sorbers formed within expanded regions is not connected 
with the DM structure (Bi & Davidsen 1997; Zhang et al. 
1998; Dave et al. 1999). Such comparison reveals also the 
potential and limitations of this approach. 

This approach relies on identification of separate ab- 
sorbers in the observed spectra what restricts the number of 
available spectra. Moreover, properties of some of the iden- 
tified lines can be distorted due to superposition of several 
lines and cannot be reliably discriminated from the influence 
of diffuse intergalactic gas (McGill 1990; Levshakov & Kegel 
1996, 1997) what introduces some additional uncertainties 
in our analysis. Our study of simulations at small redshifts 
(DDMT) shows, however, that the influence of the last effect 
depends on the power spectrum of primordial perturbations 
and for the CDM-like power spectrum, at moderate red- 
shifts, its impact is not very strong. This problem should be 
investigated in more details using the available simulations 
of absorbers. 

The redshift dependence of linear number density of 
absorbers was discussed in Demianski, Doroshkevich & Tur- 



chaninov (2000, hereafter Paper I) under the assumption 
that the neutral hydrogen traces the potential wells formed 
by DM pancakes. This model is similar, in some respects, 
to previously discussed theoretical and simulated models re- 
ferred above. In this paper we show that this model provides 
a reasonable self consistent description and interpretation 
of other observed properties of absorbers and demonstrates 
that merging of earlier formed structure elements plays an 
important role in the evolution of absorbers. This approach 
makes it possible to approximately discriminate between the 
evolution of DM structure elements and faster and randomly 
perturbed evolution of gaseous component, that is suitably 
described by the evolution of entropy of compressed gas, 
and to specify the main factors responsible for it. In par- 
ticular, we can roughly discriminate between adiabatic and 
shock formation of absorbers and show that the adiabatic 
processes may be not very important, at least for the forma- 
tion of stronger observed absorbers. Properties of observed 
absorbers demonstrate the important role of merging of ear- 
lier formed structure elements for the absorber evolution. 
Our method gives reasonable fits for the observed distribu- 
tion of Doppler parameter, b, and hydrogen column density, 
Nhi , and connects them with the basic cosmological param- 
eters, Q. m & h, and the amplitude of initial perturbations. 
Our main results are consistent with conclusions of Zhang 
ct al. (1998), Weinberg et al. (1998) and Dave et al. (1999). 
Some differences between the theoretical expectations and 
simulations will be discussed below. 

The theory cannot yet describe in details the relaxation 
of compressed matter, the disruption of structure elements 
caused by the gravitational instability of compressed DM 
and the distribution of neutral hydrogen across DM pan- 
cakes. Therefore, in this paper several parameters character- 
izing the properties of DM and neutral hydrogen distribution 
remain undetermined. They can be estimated by applying 
the discussed methods to simulations that provide an unified 
physical picture of absorber formation and evolution. 

The composition of observed absorbers is complicated 
and if at low redshifts a significant number of stronger Ly-a 
lines and metal systems is associated with galaxies (Berg- 
eron et al. 1992; Lanzetta et al. 1995; Cowie et al. 1995; 
Tytler et al. 1995; Le Brune et al. 1996) then the popula- 
tion of weaker absorbers dominates at higher redshifts and 
mainly disappears at redshift z < 2. It is not observed by 
other methods and can be associated with the population of 
weaker structure elements formed by non luminous baryonic 
and DM components and situated in extended lower density 
regions. Some number of weak Ly-a lines observed even at 
small redshifts far from galaxies (Morris et al. 1993; Stockc 
et al. 1995; Shull et al. 1996) can be considered as a trace 
of this population. 

An interesting problem arises, namely the possible re- 
construction of spatial DM distribution using the redshift 
distribution of absorbers. Two methods of such reconstruc- 
tion were proposed by Weinberg et al. (1998) and Nusser & 
Haehnelt (1998). Here we examine the method based on one 
dimensional smoothing of density field. We show that results 
depend strongly on the used sample of absorbers and on the 
method of identification of properties of DM component of 
absorbers. Interpretation of results of such reconstruction 
is now questioned and more detailed investigation of this 
problem is required. 
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Simulations of structure evolution take into account the 
impact of many important factors together and provide an 
unified picture of absorber formation and evolution. But so 
far such simulations can be performed only in small boxes 
what introduces artificial cutoffs in the power spectrum and 
makes the investigation of large scale structure evolution dif- 
ficult. Such simulations cannot yet reproduce all important 
features of interactions of small and large scale perturbations 
and the direct analysis of the observed absorbers character- 
istics might now be more perspective in this respect. Indeed, 
the large scale modulation of redshift distribution of Ly-a 
lines found by Cristiani et al. (1996) and strong nonhomo- 
geneities found at z < 2 by Williger et al. (1996), Quashnock 
et al. (1996, 1998), and Connolly et al. (1996) could be at- 
tributed to the extremely rich structure elements which are 
not yet found in simulations. 

This paper is organized as follows. The theoretical 
model of the structure evolution is discussed in Sees. 2. 
Sec. 3 contains information about the used observational 
databases. The results of statistical analysis are given in 
Sec. 4, 5 & 6. Discussion and conclusion can be found in 
Sec. 7. 



2 MODEL OF STRUCTURE EVOLUTION. 

The main observational characteristics of absorption lines 
are the redshift, z abs , the column density of neutral hydro- 
gen, Nhi, and the Doppler parameter, 6. On the other hand, 
the theoretical description of structure formation and evo- 
lution is dealing with the linear number density, n aba {z), 
temperature, T, density of DM and entropy of gaseous com- 
ponents, and with the ionization degree of hydrogen. To con- 
nect these theoretical and observed parameters a physical 
model of absorbers formation and evolution is required. 

Broad set of such models was discussed during last 
twenty years (see references above). Here we repeat some 
of the assumptions already discussed in earlier publications. 
Our consideration is based on the statistical description of 
formation and evolution of DM structure in CDM-like mod- 
els (DD99, DDMT), and it is compared with observed and 
simulated spatial distributions of DM component and galax- 
ies at small redshifts. 



2.1 Physical model of absorbers. 

Here we consider a simple self-consistent model of the ab- 
sorbers formation and evolution based on the Zel'dovich ap- 
proximation. We assume that: 

(i) The DM distribution forms an interconnected struc- 
ture of sheets (Zel'dovich pancakes) and filaments, their 
main parameters are approximately described by the 
Zel'dovich approximate theory applied to CDM-like initial 
power spectrum (DD99, DDMT). The richer DM pancakes 
can be relaxed, long-lived, and approximately stationary. 

(ii) Gas is trapped in the gravitational potential wells 
formed by the DM distribution. The gas temperature and 
observed Doppler parameter, b, trace the depth of the DM 
potential wells. 

(iii) For a given temperature the gas density within the 
wells is determined by the gas entropy created during the 



previous evolution. The gas entropy changes, mainly, due to 
the shocks heating in the course of merging of pancakes and, 
possibly, due to the bulk heating produced by local sources. 

(iv) The ionization of the gas is caused by the outer ra- 
diation field and for the majority of absorbers ionization 
equilibrium is assumed. 

(v) The evolution of observed properties of absorbers is 
mainly caused by merging, transversal compression and/or 
expansion and disruption of DM pancakes. Possible bulk 
heating of the trapped gas and possible variations of the 
intensity and spectrum of the ionizing UV radiation field can 
be also important and will provide the fine-tuning between 
the observed and expected properties of absorbers. 

(vi) In the context of the simple model we identify the ve- 
locity dispersion of DM component compressed within pan- 
cakes with the temperature of hydrogen and the Doppler 
parameter b of absorbers. We consider the possible macro- 
scopic motion within pancakes as subsonic and assume that 
they cannot essentially distort the measured Doppler pa- 
rameter. 

The formation of sheet-like DM structure elements 
(Zel'dovich pancakes) as an inevitable first step of evolution 
of small DM perturbations was certainly established both 
by theoretical considerations (Zel'dovich 1970; Shandarin & 
Zel'dovich 1989; in DD99 for CDM-like power spectrum) and 
numerically (Shandarin et al. 1995). Here we will restrict 
our consideration to the subpopulation of slowly evolving 
DM pancakes when their column density remains almost the 
same during the time comparable with H^ 1 (z). In the op- 
posite case, when rapid expansion of matter in the transver- 
sal directions takes place, the pancake is eroded and Nhi 
decreases below the observational limit. The rapid compres- 
sion transforms pancakes into filaments that is another sub- 
population of observed Ly-a absorbers. These short-lived 
pancakes can be mainly identified with a subpopulation of 
weaker absorbers with a column density Nhi < 10 13 cm~ 2 
dominated at higher redshifts z >3. 

The subpopulation of weaker absorbers also contains 
"artificial" caustics (McGill 1990) and absorbers identified 
with slowly expanding underdense regions (Bi & Davidsen 
1997; Zhang et al. 1998; Dave et al. 1999). These kinds of 
absorbers are not connected with DM structure and produce 
noise, which is stronger at higher redshifts z > 3. 

Further on, even our approximate consideration cannot 
be applied to filaments and high density clumps, when both 
the gravitational potential and the gas temperature along a 
line of sight depend essentially on the matter distribution 
across this line. So, our investigation has to be restricted 
to the subpopulation of DM composed sheet-like structure 
elements. This means that the appropriate subsample of ob- 
served absorbers has to be considered. 

Fortunately, the correlations between observed param- 
eters b and Nhi and evolutionary rate of observed linear 
number density of absorbers discussed in Paper I allows us 
to discriminate statistically the filamentary and sheet-like 
dominated subpopulations of absorbers and to perform ap- 
proximately such selection. Thus, it may be expected that 
the sheet-like DM composed absorbers dominate for b > 17 
- 20km/s and 10 13 cm" 2 < N H i < 10 14 cm~ 2 . 

Both theoretical analysis and simulations show the 
successive transformation of sheet-like elements into 
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filamentary-like ones and, at the same time, the merging of 
both sheet-like and filamentary elements into richer sheets or 
walls. Such continual transformation of structure goes on all 
the time. These processes imply the existence of complicated 
time-dependent internal structure of high density elements 
and, in particular, the essential arbitrariness in discrimina- 
tion of such elements into filaments and sheets. The mor- 
phology of structure elements can be quantitatively charac- 
terized with new powerful techniques such as the Minimal 
Spanning Tree analysis (DMRT, DDMT) and the Minkowski 
Functional (Sathyaprakash et al. 1998; Kerscher 1999). Such 
an analysis applied to observed and simulated catalogues 
demonstrates the continual distribution of morphological 
characteristics with a relatively small fraction of distinct 
high density filaments and elliptical clumps, and allows to 
estimate the degree of filamentarity and sheetness of both 
individual elements and sample as a whole. 

In this paper, as was noted above, we use the term 'pan- 
cake' for structure elements with relatively small gradient of 
properties (first of all temperature) across a line of sight. 
With such a criterion, the anisotropic halo of filaments and 
clumps can also be considered as 'pancake- like'. At the same 
time, as was discussed in Paper I, the evolutionary rate of 
observed linear number density of such absorbers can be 
similar to that typical for filaments or clumps. By imposing 
some restrictions on the observed 6 and Nhi it is possible to 
improve the statistical discrimination of these components, 
but even then the selection is not unique. More detailed 
investigation of observed and simulated properties of ab- 
sorbers is required to improve the selection criteria and to 
prepare more adequate physical model of absorbers. 



2.2 DM structure elements and Doppler 
parameter. 

Among the observed characteristics of Ly-a absorbers the 
Doppler parameter, b, is more closely linked with proper- 
ties of DM component. In this section we introduce some 
relations between characteristics of DM pancakes and the 
b parameter based both on the theoretical arguments and 
analysis of simulated DM and observed galaxy distribution 
at small redshifts. Some properties of such high density DM 
walls formed at the redshift z <Cl were analyzed in DMRT, 
DD99 and DDMT. In this Section relations between ba- 
sic characteristics of DM pancakes are introduced (without 
proofs) as a basis for further analysis. 

The fundamental characteristics of DM pancakes are 
the dimensionless Lagrangian thickness, q, and the DM col- 
umn density, fif. 



3#n 2 



(l + Zf) 8tyG 



(l + z f ) 2 q, 



(2.1) 



2.2.1 DM pancakes in Zel'dovich approximation. 

The expected probability distribution function (PDF) for 
the Lagrangian thickness of a pancake, q, can be written 
(DD99) as 



4t 2 0F 



8r 2 ' 



< £ >= 0.5 + 1/tt w 0.82, (£ 2 ) = 0.75 + 2/tt « 1.39 

where the dimensionless 'time' t(z) describes the evolution 
of perturbations in the Zel'dovich theory (Appendix A). 
More details are given in DMRT, DD99 and DDMT. 

In the Zel'dovich theory the Lagrangian thickness of 
DM pancake, q, is closely linked to the velocity of infalling 
matter, v in f, the thickness, hoM, and overdensity, 5dm, of 
compressed dark matter (DDMT): 



. . H(z)lo/3(z) 

{Vinf) « , : , q, 0~ in f 



2(1 + z) 



(M « V( h2 D M ) 



1+2 



2l T 

l + z 



y/q, 8 DM ~ -75— 



(2.3) 



2r ' 



The PDF of infalling velocity is Gaussian for Gaussian ini- 
tial perturbations (DD99, DDMT) with the mean value and 
dispersion as given by (2.3). 

Further on we will identify the kinetic energy accumu- 
lated by the DM pancake with the Doppler parameter, boM, 
and with the observed Doppler parameter, b. We will assume 
that 



°DM 



(Vinf) 



l » H2{z)q ;(qf(z) + 4T 2 [l + p(z)] 2 ). (2.4) 



12(1 + z) 2 

This assumption is valid during some time after formation 
of the pancake and allows us to reasonably describe the ob- 
served properties of absorbers. Later on, other processes 
such as the relaxation and small scale clustering of com- 
pressed matter, as well as, the pancake compression and/or 
expansion in transversal directions become important. Some 
of them will be discussed below. 

For larger redshifts, z > 2, q <JC 1, we can introduce more 
suitable notation separating out large numerical factors. We 
can take with a reasonable precision 



H{z) m H (l + zf 



267C 3/2 VQnlkm/s/Mpc, 



r w 0.06CV Z , /3«1, < = 0.25(1 + z) 



(2.5) 



what allows us to rewrite (2.3) and (2.4) more transparently 

as 



6dm « boy/e + 2£, £ = r?{l+ V 1 + V 2 ) \ V = b DM /bo, 



6.6 59.4Mpc _ no , 2 

lo « T75— h Mpc= — — — , Q m = 9Q m h, 

where f2 m is the dimensionless mean matter density of the 
universe and ffo=100 h km/s/Mpc is the Hubble constant. 
The Lagrangian thickness of a pancake, loq, is defined as an 
unperturbed distance at redshift z = between positions of 
DM particles bounding the pancake. 



(6dm) « 1.436o, {b 2 DM }~Sbl, 8 D M~y/2t, (2.6) 

6 = 33km/sC 3/2 T 2 e m 1/2 <d v , 

where £ was introduced by (2.2), and the factor 0„ ~1 de- 
scribes differences between (2.4) and (2.6). 

These relations connect the velocity dispersion within 
DM pancakes with their DM column densities and allow to 
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obtain (using the PDF (2.2) ) the expected PDF for bou as 
follows: 



N b = N q 



dq 2 e -e\/2TT 



2.2.5 Relaxation of DM pancakes. 

The analysis of simulations (DDMT) shows that in rich 
pancakes the DM particles are relaxed and gravita- 
tionally confined, and such pancakes are long-lived and 
(quasi)stationary. The relaxation of DM particles is essen- 
tially accelerated by the small scale clustering of compressed 
matter and leads to evaporation of particles with larger 
velocities, what restricts the observed Doppler parameter 
of rich pancakes. Thus, for observed and simulated galaxy 
walls, at small redshifts, b ~ 300 - 350km/s is found to be 
more typical. 

The approximate relation for the velocity dispersion of 
DM particles within relaxed pancakes can be written as fol- 
lows: 



b DM « e(z)y/{b 2 DM ) 



(2.8) 



(e) « 0.3 -0.7, 7^0.6-0.7. 



Here the random dimensionless parameter e(z) describes the 
lost of energy in the course of relaxation. Estimates of the 
factors 7 and e can be refined through a more detailed com- 
parison with simulations. 

For relaxed pancakes, instead of (2.6) & (2.7) we obtain 
that, for example, for 7 = 2/3, the expected characteristics 
of Doppler parameter can be taken as follows: 



(q 2/3 ) « 1.2(g) 2 / 3 , (b DM )*2b e, g«8r 2 (£) 



\/nb e 



exp(-£) 



g£f(vg) 

£1/6 



■ Hi) 



3/2 



3/2 



(2.9) 



Relations (2.2), (2.7) and (2.9) link the observed 
Doppler parameter with the pancake DM column density, 
q, and indicate that the distribution function of Doppler pa- 
rameter is similar to a gamma distribution. These relations 
take into account the successive merging of earlier formed 
structure elements - both filaments and pancakes - what 
leads, in particular, to formation of observed galaxy walls at 
small redshifts. 

As was discussed in Paper I, rapid expansion of pan- 
cakes in transversal directions decreases Nhi below the ob- 
servational limit of Nhi ~ 10 12 cm~ 2 . Some fraction of such 
expanded pancakes with smaller Nhi can be observed. The 
rapid compression of a pancake along one or both of the 
transversal directions decreases its surface area and also the 
probability to see such a pancake as an absorber. 

These inferences are based on theoretical arguments 
tested on simulated formation of walls at small redshifts. 
More accurate estimates of possible distortions can be ob- 
tained through comparison with representative simulations 
at high redshifts. 

As was noted above, we identify 6dm as given by (2.6) 
or (2.8) with the observed Doppler parameter, b. 



2.3 Large scale absorber distribution 

The observed spatial distribution of absorbers is very similar 
to the Poissonian distribution. Even so, the large scale mod- 
ulation of absorbers distribution is an interesting character- 
istic. If absorbers are actually linked to the DM distribution 
as was discussed in Sec. 2.2, then it can be expected that 
this modulation will be later transformed - due to the grav- 
itational instability - to the galaxy distribution observed at 
small redshifts as large and superlarge scale structure. 

The large scale matter distribution can be conveniently 
characterized by the one dimensional smoothed density field. 
The required variance of density can be expressed through 
the dimensionless moments of initial power spectrum, p(k). 
Thus, for the density smoothed over a scale r s with a Gaus- 
sian window function we have: 



(27T) 



47T 3 / 2 r s J o 



d 3 k exp[— (kr s ) 2 ] p(k) 



dk kp(k) erf (krs 



(2.10) 



where k is a comoving wave number. For r. 
1 



4vr 3/2 



dk kp(k) 



475-3/2, 



dk kp c dm(k) , 



'0 ^" • s Jo 

where p c dm is the standard CDM-like power spectrum with 
the Harrison-Zel'dovich asymptotic p{k) oc k,as k — > and 
the CDM transfer function and K c d m describes the impact 
of possible deviations of actual and CDM-like power spectra 
used for numerical estimates. For large r s we have: 

(J 2 p {r s ,z) w 18K 2 d m^ 2 (^) — , 



0.24 




(2.11) 



For such estimates we will use the measured redshift 
of absorbers and relations (2.6) and (2.9) to obtain the re- 
quired DM column density of absorbers through the mea- 
sured Doppler parameter, b. 

Such approach is similar but not identical to that used 
by Weinberg et al. (1998) and Nusser & Haehnelt (1998). 

2.4 Properties of gaseous structure elements. 

If the observed Doppler parameter characterizes the ba- 
sic properties of DM pancakes then the hydrogen column 
density characterizes the state of the gaseous component 
trapped by the DM pancakes. The gas density and the hy- 
drogen column density are sensitive to many factors. Firstly, 
the radiation field provides the high ionization and bulk 
heating of hydrogen. Secondly, the shock compression and 
heating of gas accompany the process of merging of richer 
DM pancakes. The adiabatic compression or expansion of 
gas changes also its temperature and the observed hydrogen 
column density. 

The suitable characteristic of the state of gas is its en- 
tropy. It remains constant during the adiabatic compression 
and expansion of gas and it increases due to all irreversible 
processes such as the shock and bulk heating. The entropy 
decreases only due to the radiative cooling which is usually 
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moderate. When gas temperature is rigidly bounded by the 
gravitational potential of DM distribution the entropy and 
the gas density are closely linked. 

2.4-.1 Properties of homogeneously distributed hydrogen 

The properties of compressed gas can be suitably related to 
better known parameters of homogeneously distributed gas, 
which were described in many papers (see, e.g., Ikeuchi & 
Ostrikcr 1986). In this case the baryonic density and the 
temperature can be taken as 

n b = 1.2 • Kr 5 fi(,h 2 (l + zfcm~ 3 = n ( 3 , 
n = 1.5 • KT 5 cmr 3 e 6ar , Q bar = n b h 2 /0.02, (2.12) 



T bg « 1.6 • 10 K, b bg « 16km/s, ( 



(1 + z) 



and the entropy of the gas can be characterized by the func- 
tion 

F bg = T bg /n 2 b /3 = FoC 2 , F « 2 kev • cm 2 e ta f. (2.13) 

For reference, the typical entropy of the primordial gas be- 
fore reheating, F vrm , and of gas observed in galaxies, F ga i, 
and in clusters of galaxies, F ct , are: 

F prm ~ 3 • 10~ 5 kev • cm 2 « 1.5 • 10~ 5 _F 



F ga i ~ 10 3 kev ■ cm 2 w 10 3 F 



F ci ~ 300kev • ran w 150F 

These data illustrate the range of observed variations of 
gaseous entropy. 



2.4-2 The hydrogen column density and entropy of 
absorbers. 

The observed column density of neutral hydrogen can be 
written as an integral over pancake along a line of sight 



N h 



dx 



dN ba 
dx 



-x H 



{N bar )(x H )n. 



(2.14) 



Here (N ba r) is the mean column density of baryons across 
the absorber, {xh} is the mean fraction of neutral hydrogen 
and the dimensionless parameter k characterizes the nonho- 
mogeneous distribution of neutral hydrogen across absorber. 
We will assume that both DM and gaseous components are 
compressed together and, so, the column density of baryons 
and DM component are approximately proportional to each 
other. Therefore, we can take 



(N b ar) 



n b lpq 
v{l + z) 



2.1- 10 19 cm" 2 £(&)- 



i/0„ 



{Ar){S bar ) = 0.25/oqC 1 « 0.5C(6)r z 2 C _3 Mpc. (2.15) 

Here (Ar) and {Sbar) are the mean proper thickness and 
overdensity of compressed gas above the mean density, 
£(&), TzSzC, are given by (2.2) & (2.5), v = cosip describes 
the random orientation of absorbers and the line of sight, 
and expressions (2.7) & (2.9) link functions q, b and r. 

Under the assumption of ionization equilibrium of the 



gas within a DM pancake and neglecting a possible contri- 
bution of macroscopic motions to the 6-parameter (T oc b 2 ), 
the fraction of neutral hydrogen is 



(x H ) = mC^Shar}- 



x 



(Sbar) *3 ( bbg 



3/2 



r 7 = ria • io" 12 s _1 , x = 4.1 • io" 6 e 6£ir . 

The recombination coefficient, ct rec (T), and the photoioniza- 
tion rate due to the extragalactic UV background radiation, 
T 7 , are taken as 



a rec (T) «4 - 10 



. i 'l0 4 K\ 3/4 cm 3 



, ri2~o.7 



((Black, 1981, Rauch et al. 1997). Finally, for the column 
density of neutral hydrogen we have 

6 M 3/2 (M K( 3 eH , (2.16) 



Nhi = N m (jf) 



N Q = 8.6 • 10 13 cm 



@H = 



Pi 2 2 

^bar^z 



and 5(6) was introduced by (2.2) and (2.6) or (2.9). 

The relation (2.16) links the observed column density 
of neutral hydrogen and the Doppler parameter with the 
column density of baryons and dark matter, £(&), and the 
degree of matter compression, {8 bar ). The degree of matter 
compression depends on the entropy of gas and, so, on its 
evolutionary history that allows us to discriminate between 
various models of absorbers formation and, in particular, be- 
tween adiabatic and shock compressions of the gas accumu- 
lated within absorbers. Thus, for the adiabatic compression 
of the gas which is more typical for the formation of smaller 
DM pancakes 

(Sbar) = (b/bbg) 3 , Nhi oc £(6)& 3/2 . (2.17) 

On the other hand, to describe formation of rich pancakes 
we can use the function £(&) (2.6) or (2.9). In this cases we 
have 



Nhi « 0.3iV 



Nr 



Vv 



l + ^/l + r, 2 Tia 



{Sbar} *21/4 



oljVoe -3/ 2 (5w) KC 21/4 ej/) 
Tl2 



(2.18) 



(2.19) 



for unrelaxed and relaxed absorbers, respectively. In both 
cases the expected correlation between Nhi /{Sbar) and b is 
negligible and Nhi oc {Sbar). For such pancakes the over- 
density 8dm obtained in Sec. 2.2.1 is not a good param- 
eter because the evolutionary histories and, so, degrees of 
compression of DM and gaseous components are certainly 
different. 

As seen from (2.18) & (2.19), at redshifts z > 3, C > 
1 the discussed approach can be applied mainly to rich ab- 
sorbers with Nhi > 10 13 cm~ 2 . For smaller redshifts, Q <1, 
this model can also describe the properties and evolution 
of some fraction of poor absorbers with Nhi < 10 13 cm~ 2 
connected with DM pancakes. The range of application of 
this model increases for Ti2 >1, and depends on the local 
factors which can change the parameter Qh- 

Relations (2.17) and (2.18) allow us to estimate the 
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compression factor, n{6bar)F 12 and the entropy of gas accu- 
mulated by absorbers. Two functions 

F S = (T/T bg )S^/ 3 oc b 2 (F 12 /K) 2/s , E = lnF s , (2.20) 

measure this entropy relatively to the entropy assumed for 
the homogeneous intergalactic gas (2.13). These functions 
also depend on the photoionization rate, T 7 , and the un- 
known distribution of neutral hydrogen across absorbers, 
described by the factor k. 

2.4.3 Adiabatic compression of gas. 

The velocity of infalling gas (2.3) can be rewritten more 
transparently as 

{v inf ) « 62.5km/s £r z 2 C 3/2 e- 1/2 , (2.21) 

a inf « 47fcm/s £ 1/2 r z 2 C 3/2 e- 1/2 , (2.22) 

and for poor pancakes with £ < 0.25£ 3,/2 , when (vi„f) < 
bt g ~ 16km/ s, the formation of DM pancakes is accompa- 
nied by the adiabatic compression of gas. As seen from Eq. 
(2.17), in this case 

£(&) oc N HI b-' i/2 < 1, 

and Eq. (2.2) shows that the homogeneous distribution can 
be expected for both £ and N H ib~ 3 ^ 2 . 

This conclusion can be, in principle, tested using the 
observed sample of absorbers. Unfortunately, the available 
sample of weaker Ly-a lines is poor and incomplete and can- 
not provide the statistics required for such test. Moreover, 
such absorbers are more strongly influenced by random lo- 
cal factors that also can destroy the expected correlation 
between b and Nhi- On the other hand, the same factors 
can generate the local shock waves without any connection 
with the evolution of dark matter. Moreover, this sample 
is complex and it also includes absorbers which are not a 
byproduct of the process of DM structure formation. 

2.4-4 Shock compression of gas. 

Under the shock compression the gas density increases not 
more than 4 times while the Doppler parameter increases 
proportionally to v in f . Simple estimates similar to that given 
in Zel'dovich & Novikov (1983) for "Zel'dovich pancakes" 
and Meiksin (1994) for the collapse of ionized gas into slabs 
show, however, that even for the strong compression of ho- 
mogeneously distributed gas the smooth profile of Vi n f oc £ 
(2.21) leads to strong adiabatic compression of the gas be- 
fore formation of shock waves (see, e.g., Nath 1997). 

For the observed range of Doppler parameters the most 
plausible scenario of gas evolution is shock compression of 
the gas already accumulated within the DM confined pan- 
cakes and clouds ('grain' model). This process is typical for 
the merging of earlier formed structure elements. Such mat- 
ter distribution reduces the adiabatic compression of the 
gas during the stage when clouds approach each other. The 
shock wave is formed at the moment of merging of clouds and 
only shock compression occurs. This process results in an es- 
sential increase of the entropy because of the limited growth 
of density and larger growth of the temperature T oc v 2 n f . 



This scenario can be realized if a high matter concen- 
tration within DM confined structure elements has occurred 
before creation of observed absorbers. Such strong matter 
concentration within clouds, similar in some respects to that 
considered in the 'mini-halo' model (Rees 1986, Miralda- 
Escude & Rees 1993), is consistent with the theoretical ex- 
pectations (DD99) because, for the CDM-like transfer func- 
tion and Harrison - Zel'dovich primordial power spectrum, 
the majority of matter should be accumulated by low mass 
clouds with Mdm ~ 10 7 MQ(£l m h 2 )~ 2 already at redshifts 
z > 5. This conclusion is consistent with simulations (see, 
e.g., Zhang et al. 1998) which demonstrate that only ~ 5% 
of baryons remains in a smoothly distributed component. 

The analysis of observational data (see below) shows 
that for observed samples there is a strong correlation be- 
tween the entropy and the Doppler parameter, while b and 
Nhi are practically not correlated. For richer absorbers, 
with Nhi > 10 13 cm -2 , the entropy distribution can be ap- 
proximated as Fs oc b 2 what agrees well with this model of 
shock heating of gas with strongly nonhomogeneous distri- 
bution. 

The entropy, E, given by (2.20) is an additive func- 
tion which accumulates the successive contributions of shock 
and bulk heating during all evolutionary history of a given 
gaseous element. If the shock heating of gas dominates, then 
the growth of entropy, S, can be described as a random pro- 
cess - similar to the Brownian motion - with a successive 
uncorrelated jumps of entropy at each step. This means that 
the expected PDF of the entropy, E, should be similar to 
Gaussian. 

For the samples of poorer absorbers with Nhi < 
10 13 cm -2 there are both substantial correlations between 
b and Nhi and between the entropy and the Doppler pa- 
rameter. This fact indicates the complex character of such 
absorbers. If they can be associated with pancake-like DM 
elements then for some fraction of such absorbers the adia- 
batic compression or bulk heating could be significant, while 
some of them could be formed due to expansion of richer ear- 
lier formed pancakes. Some fraction of such absorbers can be 
also related to artificial pancakes discussed by McGill (1990) 
and Levshakov & Kegel (1996, 1997) and to absorbers situ- 
ated within "minivoids" what increases uncertainty in their 
discrimination. 

2.4-5 Bulk heating and cooling of the gas. 

The bulk heating and cooling of the compressed gas can be 
also essential for the evolution of properties of absorbers. 
The main factor is the random spatial variation of the 
spectrum of ionizing UV radiation field generated by local 
sources (Zuo 1992, Fardal & Shull 1993). The action of this 
factor can be characterized by the mean energy injected at 
photoionization, T 7 w(5 - 10) -W 4 K for suitable spectra 
of UV radiation (Black 1981). In extremal cases such vari- 
ations are observed as a proximity effect (Bajtlik, Duncan 
& Ostriker 1988). In the model of DM confined absorbers 
discussed here, when the gas temperature is given by the 
gravitational potential of DM component, the bulk heating 
changes the density and entropy of compressed gas. 

The influence of the bulk heating can be enhanced by 
the pancake disruption due to the clustering of both DM and 
baryonic components. This process is usually accompanied 
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by adiabatic reduction of temperature in expanded regions 
when the role of the bulk heating becomes more important. 
As is shown in Appendix B, in the general case, for redshift 
dependent T(z) & T 7 (z) we have 



p 3/2 



p 3/2 



(*/) + 



Cts I 

J z 



dx 



T(x) 



H(x) (l + x)T bg 



a s — 1.5h £ Qbar, 



(2.23) 



where the entropy function Fs was introduced by (2.20) and 
Zf is the redshift of pancake formation. For example, un- 
der condition of slow variation of the DM distribution when 
T(z) w const, and for T 7 w const, we have 



F 3 s /2 (z) = F^(z f ) + a s H [t(z) 



r 3/2 



Ibg 



t(z) 



f 

•J z 



dx 



H(x)(l + x) 



that indicates the slow growth of entropy and drop of the 
density of absorbers. 

At higher redshifts the proper separation of structure 
elements is about D sep ~0.5 - l/i _1 Mpc that is compara- 
ble with sizes of galactic halos observed at small redshifts 
~ 0.2/i _1 Mpc (see, e.g., Bahcall et al. 1996). This means 
that the evolution of structure elements at such redshifts 
could be interdependent and activities in galaxies can in- 
crease the gas entropy within nearby pancakes. Analysis of 
Cen & Ostriker (1993) shows that the input of explosive en- 
ergy can be essential in the vicinity of virialized objects. In 
this case the gas overdensity decreases in proportion to the 
injected energy and the additional correlations between the 
Dopplcr parameter and column density of neutral hydrogen 
do not appear. The observations of metal systems together 
with relatively weak Ly-a lines shows that sometimes the 
influence of such heating can be important. 



2.4-6 Evolution of hydrogen column density 

The observed column density of absorbers is changed due 
to irreversible processes such as shocks and bulk heating, 
and due to adiabatic expansion and/or compression of pan- 
cakes caused by the transversal motions of DM component 
(DD99, Paper I). The heating of the gas does not distort the 
DM distribution and the ionization equilibrium and, as was 
discussed in Sec. 2.4.5, also the relatively slow decrease of 
Nhi oc (S 2 ar ). The adiabatic compression and expansion of 
DM and gaseous components also does not distort the ion- 
ization equilibrium, but changes the observed column den- 
sity Nhi oc 6 9 ^ 2 , and can lead to the appearance of rare 
absorbers with larger and smaller Nhi and b. 

These simple arguments demonstrate the important 
role of adiabatic evolution of DM pancakes for the formation 
of observed sample of absorbers. Together with the merging 
these processes lead to the fast evolution of linear density of 
absorbers discussed in Paper I and are essential for the for- 
mation of both rare absorbers, with larger Nhi, and weaker 
absorbers. 



2.5 Theoretical expectations 

The previous consideration can be shortly summarized as 
follows: 

(i) For the DM confined absorbers the observed Doppler 
parameter, b, is closely linked to the column density of DM 
pancakes, q, and the expected PDFs of the b parameter, 
(2.7) and (2.9), are similar to the gamma-distribution. The 
density field can be described as a set of discrete pancakes 
with masses oc q(b). 

(ii) For such absorbers only weak correlation between the 
Doppler parameter, b, and the column density of neutral hy- 
drogen, Nhi, is expected under the assumption of ionization 
equilibrium. 

(iii) The same factors result in the strong correlation be- 
tween the entropy of compressed gas and the Doppler pa- 
rameter. Due to merging and shock heating the gas entropy 
increases. This growth can be described as a random process 
with the function E = In Fs increasing discontinuously with 
jumps depending on the mass of merging pancakes. The pos- 
sible bulk heating caused by the local sources results also in 
random jumps of entropy. More detailed description of en- 
tropy evolution can be obtained on the basis of Fokker-Plank 
equation that implies however more detailed description of 
the pancake evolution and assumptions about the local ac- 
tivity of galaxies. In any case approximately Gaussian dis- 
tribution of £ can be expected. 

(iv) As is seen from equations (2.18), (2.19) & (2.20) for 
such pancakes the approximately Gaussian PDFs are also 
expected for log Si, ar and log Nhi- 

(v) The adiabatic compression of the gas dominates for 
poorer DM pancakes with Nhi < W 13 cm~ 2 and creates 
much stronger correlation between Nhi and b. This corre- 
lation can be, however, partly destroyed by the action of 
local random factors. Formation of such absorbers is often 
not accompanied by formation of DM structure elements. 

As is seen from (2.16) and (2.20) the observed estimates 
of entropy, £, and the column density, Nhi, depend on the 
intensity of UV ionizing radiation (factor Ti2). This means 
that local variations of the ionizing rate introduce the essen- 
tial random factor in observational estimates of both entropy 
and column density. The rapid growth of observed number 
of galaxies at z < 3 (Steidel et al. 1998; Giavalisco et al. 
1998) correlates well with the observed rapid evolution of 
linear density of absorbers, n a bs, (Paper I) what points to- 
ward galaxies as the essential factor of observed evolution of 
absorbers. 



2.6 Characteristics of absorbers in simulations 

High resolution simulations of dynamical and thermal evo- 
lution of gaseous component cited above give us valuable 
examples of absorbers formed at redshifts z <5 and allow to 
clarify many peculiarities and properties of these processes. 
In particular, they confirm the existence of high density 
clumps, filamentary and sheet-like components of structure 
at high redshifts and demonstrate that these elements accu- 
mulate up to 95% of baryons and reproduce well the main 
observed characteristics of absorbers. They confirm the dom- 
ination of adiabatic compression of baryonic component in 
less massive pancakes and show that some fraction of weaker 
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absorbers can be identified with underdense structures (Bi 
& Davidsen 1997; Zhang et al. 1998; Dave et al. 1999). List 
of these examples can be essentially extended. 

At the same time, the abilities of simulations in statisti- 
cal description of absorbers evolution are limited as the main 
simulations are performed with small box sizes (~10 Mpc) 
that is comparable with the mean separation of galaxy fila- 
ments at small redshifts, and is smaller then the mean sep- 
aration of richer galaxy walls (~40 - 60/i -1 Mpc) and even 
their Lagrangian size (~15 - 20ft~ 1 Mpc). The small box 
sizes restrict the important influence of large scale pertur- 
bations on the evolution of small scale structure and does not 
allow to establish direct connection of structure at high red- 
shifts with galaxy distribution observed at small redshifts. 

Moreover, some results obtained in simulations should 
be explained in more details. Thus, for example, in all papers 
the very important problem of quantitative characteristics 
of absorbers morphology is not discussed and the merging 
of pancakes and filaments is not considered. Such merging is 
certainly responsible for the formation of observed galactic 
walls and was discussed in DD99, DDMT, Paper I and in 
Sec. 2.2 & 2.3 above as one of the main factors of absorbers 
evolution. Its action can be easily traced using the analy- 
sis of entropy of compressed gas. The distribution function 
of temperature (Fig. 10 in Zhang et al. 1998) is quite dif- 
ferent from the observed distribution of Doppler parameter. 
The strong contribution of macroscopic motions in simu- 
lated Doppler parameter plotted in Fig. 13 increases it about 
of 2 - 5 times with respect to the thermal value what prob- 
ably implies the supersonic character of typical macroscopic 
motions within simulated structure elements. In contrast, 
the analysis of DM evolution in large boxes demonstrates an 
approximate isotropy of subsonic velocity dispersion within 
DM structure elements at small redshifts, and its close con- 
nection with the measured richness of such element (DMRT, 
DDMT). The physical reasons of the power distribution of 
Nhi plotted in Fig. 15 of the same paper, in the range of 5 
orders of magnitude, are also not explained. 

Even these examples demonstrate an essential scatter of 
some quantitative characteristics of structure obtained with 
different codes and parameters of simulations (see also dis- 
cussion in Melotte et al. 1997; Splinter et al. 1998; Theuns 
et al. 1999). They show also the usefulness of more detailed 
comparisons and tests of consistency of results obtained with 
different approaches and box sizes, and based on larger set 
of characteristics. Bearing in mind the relatively small num- 
ber of such complicated simulations, these differences, list of 
which can also be continued, seem to be natural. 



3 THE DATABASE. 

The present analysis is based on the spectra available in the 
literature. The list of the used sources of data is given in 
Table 1. As was discussed in Paper I absorbers with 6 > 
17km/s or \og(Nni) < 14 can be probably related to the 
sheet-like component of structure. In spite of statistical char- 
acter of such discrimination it allows us to obtain more ho- 
mogeneous sample of sheet-like absorbers used for compari- 
son with theoretical expectations. The list of available Ly-a 
lines was arranged into three samples. Two of them, Q12 and 
Q14 include 2073 and 2378 lines from the first 12 and all 14 



Table 1. QSO spectra from the literature 



Name 


2em 




Zmax 


FWHM 


No 










km/s 


of lines 


0000 - 260 1 


4.11 


3.4 


4.1 


7 


431 


0014 + 813 2 


3.41 


2.7 


3.2 


8 


262 


0956 + 122 2 


3.30 


2.6 


3.1 


8 


256 


0302 - 003 2 


3.29 


2.6 


3.1 


8 


266 


UDoD -+" DoU 


Q 17 
O.l/ 


Z.O 


o.U 


Q 
O 


olo 


1946 + 766 3 


3.02 


2.4 


3.0 


8 


461 


0055 - 259 4 


3.66 


2.9 


3.1 


14 


313 


2126 - 158 5 


3.26 


2.9 


3.2 


11 


130 


1700 + 642 6 


2.72 


2.1 


2.7 


15 


85 


1225 + 317 7 


2.20 


1.7 


2.2 


18 


159 


1101 - 264 s 


2.15 


1.8 


2.1 


9 


84 


1331 + 170 9 


2.10 


1.7 


2.1 


18 


69 


1033 - 033 10 


4.50 


3.7 


4.4 


18 


299 


2206 - 199 11 


2.56 


2.1 


2.6 


11 


101 



1. Lu et al. (1996), 2. Hu et al., (1995), 3. Kirkman & Tytler 
(1997), 4. Cristiani et al. (1995), 5. Giallongo et al. (1993), 6. 
Rodriguez et al. (1995), 7. Khare et al. (1997), 8. Carswell et al. 
(1991), 9. Kulkarni et al. (1996), 10. Williger et al. (1994), 11. 
Rauch et al. (1993), 

spectra, respectively, under conditions that b > 17km/s and 
10 13 cm~ 2 < Nhi < 10 14 cm~ 2 . First condition allows us to 
discriminate and to exclude absorbers which could be pos- 
sibly formed due to the adiabatic compression. 

The sample Q0612 contains 469 lines with Nhi < 
10 13 cm~ 2 and b < 30km/s from the first 6 QSOs. This sam- 
ple is incomplete as was discussed by Hu et al. (1995). It can 
be used for the estimates of properties of poor absorbers. 

The line distribution over redshift is nonhomogeneous 
and the majority of lines are concentrated at 2 w 3. The 
distribution of absorbers at z > 3.2 is based mainly on the 
spectrum of QSO 0000-260 (Lu et al. 1996) and here the 
line statistics is insufficient. The inclusion of the spectrum 
of QSO 1033-033 extends the redshift interval up to z « 4.4, 
but cannot eliminate small representativity of the sample at 
such redshifts. 



4 STATISTICAL CHARACTERISTICS OF 
ABSORBERS 

Some statistical characteristics of absorbers were found for 
samples Q12 and Q14 for four ranges of redshifts. Main re- 
sults are listed in Table 2 and plotted in Figs. 1-4. 



4.1 Distribution of observed Doppler parameter, 
b, and DM column density, q 

As was discussed in Sec. 2.2, the observed distribution of 
Doppler parameter, b, can be interpreted as the distribution 
of DM column density of absorbers. To check this hypoth- 
esis the observed distribution of b/{b) was fitted to the two 
parameters function similar to (2.6) & (2.7) 

N b = Cl erf(,/y)^^ exp(-y), (4.1) 



© 0000 RAS, MNRAS 000, 000-000 



10 Demianski & Doroshkevich 



Table 2. Main parameters of Ly-a absorbers. 







Zmax 




c 


(b) 
km/s 


(9> 


C2 


T 6 


r bH 


(Fs) 


Pbs 
















Ql2 














a 


1.7 


4.1 


2073 


0.98 


36.9 


0.2 


1.6 


0.17 


0.04 


4.8 


1.89 


0.9 


b 


2.0 


3.0 


1118 


0.92 


37.6 


0.2 


1.5 


0.17 


-0.04 


5.0 


2.14 


0.9 


c 


2.5 


3.5 


1392 


0.98 


36.8 


0.2 


1.9 


0.18 


0.02 


4.9 


2.05 


0.9 


d 


3.0 


4.1 


813 


1.10 


36.0 


0.2 


2.0 


0.19 


0.16 


4.8 


1.59 


0.8 




1.7 
2.0 
2.5 
3.0 


4.4 
3.0 
3.5 
4.4 


2378 
1190 
1403 
1046 


1.00 
0.92 
0.98 
1.11 


38.2 
37.3 
36.8 
39.5 


Ql4 
0.2 
0.2 
0.2 
0.3 


1.9 
1.6 
1.7 
1.8 


0.16 
0.17 
0.17 
0.21 


0.11 
-0.02 
-0.02 

0.23 


6.1 
4.8 
4.6 
5.4 


1.91 
2.10 
2.05 
1.86 


0.8 
0.8 
0.8 
0.8 



The correlation coefficient r^u, exponent pt, a 
4.3 and 4.2, respectively. 



and the entropy dispersion erg are introduced in Sees. 




Figure 1. The distributions of Dopplcr parameter, b/(b), for the 
sample Q12 of Ly-a lines for different redshifts. The best fits (4.1) 
are plotted by dashed lines. The main parameters of subsamplcs 
are listed in Table 2. 




0.0 



0.5 



1.5 



2.0 



1.0 

q/<q> 

Figure 2. The distributions of DM column density, q/ (q), for the 
sample Q12 of Ly-a lines for different redshifts. The best fits (2.2) 
are plotted by dashed lines. The main parameters of subsamplcs 
arc listed in Table 2. 



y 



\b 2 /(bf 



1 + y/l + 



for both samples and four different ranges of redshifts. For 
each absorber the dimensionless column density of DM com- 
ponent, q, was found using Eq. (2.6) and its distribution was 
fitted to a two parametric function 

N q « c 3 e-" , y = c A q/(q) . (4.2) 

The main results are listed in Table 2 and plotted in Figs. 1 
& 2. 

The observed properties of the Doppler parameter, 6, 



are found to be surprisingly stable and independent of sam- 
ples or redshift ranges. In all cases we have 



(b) = 37km/s, o b = 0.55(6). 



(4.3) 



The functions (4.1) & (4.2) fit well the observed blkq 
distributions for all samples under consideration. Fit pa- 
rameters C1&C2 and C3&C4 describe the cutoff in observed 
PDFs at 6 w 0.5(6). This cutoff naturally appears for the 
Doppler parameter of gaseous component as in both sam- 
ples, Q12 & Q14, absorbers with 6 < 17km/s m 0.5(6) were 
excluded from the analysis. This cutoff increases the ob- 
served (6) and (q) by the factor of C2 and C4, respectively, 
in comparison with theoretical expectations for the velocity 
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-2 



-1 



Figure 3. The distribution function of absorbers 'entropy', S, 
for the sample Q12 for four different redshifts. The best Gaussian 
fits are shown by dashed lines. 



and column density of DM component. Some deficit of ab- 
sorbers with larger 6 > (2 - 2.5)(6) is seen in Fig. 1. For the 
PDF (4.1) the expected value is a b = 0.7(6). The difference 
between observed and expected a b is also explained by the 
same cutoff and for y > 0.5, the expected a b — 0.6(6} is 
practically identical to (4.3). 

The fitting parameters C2&C4 allow us to correct mea- 
sured (6) and (q) for the cutoff at 6 w 0.5(6} and to estimate 
the amplitude and time scale of the DM structure evolution, 
t, as given by (2.2) and (2.6): 



overdensity (S bar ), of gaseous halo because they depend on 
the unknown parameter k introduced in (2.14) which charac- 
terizes the distribution of neutral hydrogen across absorber. 

4.2 Entropy of absorbers. 

The observed characteristics of the entropy described by the 
functions Fs and E (2.20) were found for the same subsam- 
ples of absorbers and the main results are plotted in Figs. 3 
and listed in Table 2. The function Fs can be fitted to the 
expression 

Fs = (Fs) exp(S), S = X - p b Jn ( , (4.6) 



(b) 



n 



{\ + z)t = (l + z) 



(g) 
6.55C4 



which discriminates between the regular variations of the 
entropy described by the first terms in (4.6), and the inte- 
gral action of random factors described by the function S. 
The PDF of S plotted in Figs 3 is fitted to Gaussian func- 
tions with as listed in Table 2. It is found to be sufficiently 
stable and weakly sample depended. At larger S the deficit 
of absorbers with higher entropy is seen as an asymmetric 
shape of PDF. It can be caused by the deficit of observed 
absorbers with larger 6 and smaller Nhi- Estimates of (Fs) 
listed in Table 2 show that the entropy of absorbers with 
larger 6 > (4 — 5) (6) is similar to the entropy of gas accu- 
mulated by clusters of galaxies as is given in Sec. 2.4.1. 

The value p bs ~2 shows that merging of pancakes is the 
main factor of entropy evolution. It shows also that among 
pancake merging there is one which provides the main jump 
in both the Doppler parameter and entropy and their cor- 
relation. The impact of smaller jumps in the course of the 
merging of pancakes, random local variations of ionizing UV 
radiation, a possible bulk heating and other random factors 
is well described by the Gaussian distribution of S-function. 
It agrees well with negligible correlation of measured Nhi 
and 6, since the correlation coefficient r b n <0.1. 

The small range of observed redshifts, 0.8< £ <1.2, does 
not allow more detailed analysis of possible variations of the 
mean entropy with z. Nonetheless, the weak decrease of p bs 
and corresponding growth of r b H at 2 > 3 shows that for 
larger z the merging of pancakes could be accompanied by 
other irreversible processes and/or by adiabatic gas com- 
pression. 



= (1+*: 



(b) 



720km /s c 2 



C 3/4 6i /4 w (0.17 ±0.1). (4.4) 



C2 and r b listed in Table 2 are found to be weakly depen- 
dent on redshift z and on subsamples used. This value r b is 
roughly consistent with the values (1 + z)r ~0.2 - 0.4 ex- 
pected for low density cosmological models (Appendix A). 

We can also estimate the typical proper thickness of 
DM absorbers, hoM, which is linked with r by Eq. (2.3), we 
have 

lu \ 8 , t2 130k P c ( n V fA « 

<fc DM )«_J -_~— — — . (4.5) 



This value is approximately consistent with sizes found by 
Dinshaw et al. (1995) and is similar to the measured sizes 
of galactic halos at small redshifts (Lanzetta et al. 1995; 
Bahcall et al. 1996) and to results obtained in Paper I. We 
cannot obtain reasonable estimates of the size, (Ar), and 



4.3 Observed distribution of column density of 
HI. 

One of the most enigmatic feature of the distribution of ob- 
served column density of neutral hydrogen is that it can 
be approximated by a single power law with a power index 
/3 H = 1.5 ± 0.05 in the range 10 13 cm" 2 < N H i < 10 22 cm~ 2 
(see, e.g., Tytler 1987, Hu et al. 1995, Kim et al. 1997). 
For smaller Nhi the deviations from the fit are usually as- 
signed to the incompleteness of samples. It is interesting that 
similar power distribution was found (but not explained) in 
simulations as well (Zhang et al. 1997, 1998) 

In the considered model of absorbers, when the gas tem- 
perature and the Doppler parameter, 6, are hardly linked 
with characteristics of DM distribution, the well known neg- 
ligible correlation of Nhi and 6 indicates the weak connec- 
tion of Nhi with the properties of DM component of ab- 
sorbers, such as, in particular, the DM surface density of 
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Figure 4. The distribution function of log Nhi for the sample 
Ql2 for four ranges of rcdshifts. Parameters of samples and fitting 
parameters are listed in Table 2. The best fits (4.8) and (4.9) are 
plotted by dashed and dotted lines, respectively. 



pancake, q. For samples under consideration the correlation 
coefficient 



{logNni * log b) — {log N hi) {log b) 

OlogN HI Olog b 



(4.7) 



was estimated directly using measured b and logNni and in 
all the cases r bH <0.1 were found (Table 2). For the DM 
dominated absorbers the column densities of baryonic com- 
ponent and, therefore, also Nhi depend mainly on the en- 
tropy of gas, Fs or E, which accumulates the contribution 
of irreversible processes during all evolutionary history of 
the compressed gas. The action of this factor disconnects 
the surface densities of baryonic and DM components, Nhi 
and qlkb. Essential variations of the entropy implies essential 
variations of observed Nhi even for the same temperature, 
6, and DM surface density, q. This means that it is difficult, 
if not impossible, to explain this joint power distribution in 
such a wide range of column density. 

The negligible correlation between the observed 6 and 
Nhi implies also that 



Nhi oc {S bar ) oc {F s /V) 



2,-2/3 



log Nhi 



-S + const. 



This means that at least for the considered subsample of 
pancake- like absorbers the PDF for observed Nhi should be 
similar to PDF of the function S, discussed in previous Sub- 
sec. Therefore, the expected PDF of log{Nni /N m ), N m = 
10 13 cm -2 , is Gaussian and can be written as 



Nh = No exp 



2* H 



loi 



.2 ( Nhi 
' V N 



(4.8) 



The distribution function of log(iVijj) is plotted in Fig. 
4 for four redshift intervals, for the sample of 12 QSO and 
for all observed lines, with Nhi > 10 13 cm~ 2 . The physical 
model of absorbers evolution introduced in Sees. 2.2 and 
2.4 applies only to absorbers with Nhi < 10 14 cm~ 2 . The 
distribution of stronger absorbers, which are probably linked 
with filaments and high density clumps, should be discussed 
in a context of spatial distribution and evolution of these 
components of structure. The statistics of absorbers with 



Nhi > 10 1 



is limited, and the right part of Figs. 4 



is rather illustrative but, even so, these Figs, demonstrate 
that the relation (4.8) also approximately fits the observed 
distribution of Nhi up to Nhi ~ 10 17 cm~ 2 . 

The power distribution of Nhi corresponds to the ex- 
ponential distribution of log{Nni) and can be written as 



Nh oc exp 



2.30(1 



Uog( 



N, 



N n 



(4.9) 



This fit is plotted in Figs. 4 by dotted lines for /3h ~ 1.5 (Hu 
et al. 1995; Kim et al. 1997). It agrees with the observed PDF 
Nh for larger log{Nni) but predicts some excess of weaker 
lines. 

These results show that the problem deserves further 
investigation in a wider range of redshifts with a more rep- 
resentative sample of absorption lines. 



5 PROPERTIES OF WEAKER ABSORBERS 

The sample of observed weaker absorbers with Nhi < 
10 13 cm -2 and b < 30km/s is incomplete and composed of 
only 469 lines in 6 QSOs listed in Table 1. This factor as well 
as the much more complicated composition of this subpopu- 
lation indicate illustrative character of our consideration in 
this subsection. 

As was discussed in Sec. 2.2 such absorbers can be cre- 
ated by adiabatic compression of homogeneous gas or the 
disruption and adiabatic expansion of earlier formed rich 
pancakes. For the first subpopulation the entropy is not 
changed and « 1 can be expected. For such pancakes 
b and Fs are evidently uncorrelated, but significant corre- 
lation between Nhi and b can be expected. For the second 
subpopulation Fs > 1 and, as before, it correlates with b, but 
Nhi and b are only weakly correlated. Moreover, for such 
absorbers the contribution of artificial pancakes discussed in 
McGill (1990) and Levshakov & Kegel (1996, 1997) and ab- 
sorbers formed within "minivoids" (Zhang et al. 1998) can 
be more significant. 

This means that the population of weaker absorbers 
consists of objects with different evolutionary history, and 
the available observational data do not allow us to discrim- 
inate between these subpopulations. So, we have to restrict 
our consideration to the analysis of PDF for the column den- 
sity of DM component q oc Nnib^ 3 ^ 2 , as was described in 
Sec. 2.4.3. 

This PDF, N q , is plotted in Fig. 5. As was expected for 
smaller q, N q « const, especially at larger redshifts z > 3. 
The significant correlation between log Nhi and log b, since 
rbH ~ 0.4, agrees also with the probable substantial contri- 
bution of adiabatically compressed subpopulation of weaker 
absorbers. The 'tail' of absorbers with larger q/{q) can be 
associated with the subpopulation of disrupted and/or ex- 



© 0000 RAS, MNRAS 000, 000-000 



Properties of observed Ly-a forest 13 




0.6 01 

q/<q> 

Figure 5. The distributions function of q oc Nnib^ 3 / 2 for the 
sample of weaker absorbers with Njji < 10 13 cm~ 2 and b <30 
km/s for the same ranges of redshifts. 



panded earlier formed richer absorbers. The fraction of such 
absorbers increases for smaller z. 



6 SPATIAL DISTRIBUTION OF ABSORBERS 

The observed redshift distribution of absorbers contains 
significant information about the spatial matter distribu- 
tion at high redshifts, as suggested in particular by Oort 
(1981, 1984). The available absorption spectra cover typi- 
cally the range D ~ 200 - 300(fi m )~ 1/2 /i _1 Mpc, whereas a 
mean separation of weak absorbers with Nhi ~ 10 12 cmT 2 
is ~l(r2 m )~ 1/ ' 2 /i~ 1 Mpc. This allows us to analyze both the 
small and large scale matter distribution. 



6.1 Small scale absorber distribution 

During last years weak clustering of Ly-a absorbers on small 
scales (Av < 300km/s) have been found for a few quasars 
(Webb 1987; Cristiani et al. 1995; Hu et al. 1995; Ulmer 
1996; Fernandes-Soto 1996). For many other objects such 
clustering is negligible and the absorbers distribution is 
nearly Poissonian. Stronger small scale correlation, found 
for metal lines, is naturally explained since several lines can 
be generated in the same gaseous cloud. At the same time 
much stronger correlation of galaxies is found both at small 
and even high redshifts (see, e.g., discussion in Steidel et al. 
1998; Governato et al. 1998; Giavalisco et al. 1998)). 

This divergence can be naturally explained in the frame- 
work of discussed here approach, when absorbers are asso- 
ciated with structure elements - filaments and pancakes - 



rather then with galaxies. The distribution of structure ele- 
ments along a random straight line is expected to be Poisso- 
nian - like (DD99, DDMT) what is consistent with the ob- 
served distribution of filaments and walls at small redshifts 
(Doroshkevich et al. 1996). In contrast, the point - like high 
density clumps, which can be associated with 'galaxies', arc 
mainly incorporated into filaments and massive pancakes. 
Even for randomly distributed filaments and pancakes this 
concentration introduces some regularity in the spatial dis- 
tribution of such clumps in comparison with the 3D Poisso- 
nian distribution. As was shown by van de Weygaert (1991) 
and Buryak & Doroshkevich (1996) such point concentra- 
tion generates the 3D correlation function similar to that 
observed for galaxies. 

Similar situation occurs, for example, for the LCRS, 
where the usual correlation function of galaxies was found 
(Tucker et al. 1997). However, for the same observed sample 
the ID distribution of both filaments and walls (in radial 
and transversal directions) were found to be Poissonian - 
like. 



6.2 Large scale modulation of absorbers 
distribution 

The smoothed absorbers distribution can be compared with 
expectations discussed in Sec. 2.3. Three expressions for 
DM distribution are examined: 1) as was discussed in Sec. 
2.2.2, for relaxed pancakes q oc 6 3//2 can be expected, 2) 
if the compre ssed matter is not completely relaxed, then 
goc-^/l + j) 2 — lis more probable, and 3) for comparison 
with the previous definition, the hypothesis q oc Nhi is also 
tested. The third expression is correct for weaker absorbers 
discussed in Sec. 2.4.3 and Sec. 5 and, possibly, can be also 
applied to high density peaks associated with galaxies but 
its application to main part of absorbers is in question. 

For separate QSO the observed absorbers distribution 
along a line of sight was averaged with the Gaussian window 
function and the smoothed density was found as 



<P> W n (x) ^ (q) 



Vi), 



W{y-x) : 



27rr, 



exp 



2r 2 



f 



(6.1) 



dyW(x-y), 



where ythqt are coordinates and DM column density of ab- 
sorber, and W n is a normalization factor taking into account 
the finite size of observed spectra. The results obtained for 
absorbers with Nhi > 10 13 cm~ 2 and averaged over 15 r s , 
7.5fi m 1/2 /t- 1 Mpc < r s < 40n m 1/2 ft" 1 Mpc, are plotted in 
Fig. 6 together with the dispersions plotted as error bars. 
Averaging over 12 QSO gives for the three used definitions 
of q: 

<(l + z)r(z))^(0.11±0.06)(e m ) 1/4 , 



1/4 



((1 + z)t(z)) « (0.12±0.05)(6 m ) 
((1 + z)t(z)) » (0.40 ± 0.22)(9 m ) 1/4 . 



(6.2) 



where, as before, O m = 9Q. m h 2 . The large scatter of points 
shows that the used statistics of absorbers is insufficient and 
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Figure 6. Redshift dependence of the smoothed amplitude, 
(1 + z)t, vs. z for the sample of 12 QSO. Three definitions of DM 
density are used: pr,M °t h 3 / 2 (top panel), Pdm °^ (1 + t? 2 ) 1 / 2 — 1 
(middle panel), and pum °^ Nhi (bottom panel). Open squares 
show the same amplitudes obtained for a random redshift distri- 
bution of the same absorbers. The mean values are plotted by the 
dashed lines and the dispersions are shown by dotted lines. 

does not allow to obtain reasonable estimates of the ampli- 
tude. 

These results show that first and second definitions of 
DM column density lead to similar estimates of t, whereas 
the third definition increases it about 4 times. 

It is interesting to compare these results with similar 
estimates obtained for randomly distributed absorbers with 
the same b and Nhi- Results averaged over 100 random 
realizations are plotted in Fig. 6. For the first and second 
definitions of DM column density, q, relative to the mean 
value, does not differ by more than a factor of 2 - 3 and 
the use of random redshift decreases the dispersion about 
2 times. For the third definition variations of q around the 
mean value are much larger and can exceed a factor of ~ 
100. Such a large discrepancy is a result of small fraction of 
absorbers with largest Nhi- 

This estimates of r q £s Tb differ by about 2-3 times 
in comparison with the values r b obtained in Sec. 4.1. This 
divergence can be caused by various factors, the most impor- 
tant are probably the small representativity and relatively 
small size of spectra used, what increases the uncertainty of 
all estimates for larger smoothing scales. 

These results, as well as estimates of Tb presented in 
Sec. 4.1, demonstrate the limited abilities of investigations 
of properties of absorbers with respect to the description 
of DM distribution. Both theoretical analysis and investi- 
gations of simulations (Governato et al. 1998; Jenkins et al. 
1998) show that at redshifts z ~ 3 high density filaments ac- 
cumulate main fraction of galaxies and are the most typical 
and conspicuous structure elements. But such filaments as 
well as separate galaxies constitute relatively small fraction 
of the observed absorbers because of their relatively small 



size, and are represented mainly by the lower density halo. 
This means that the observed absorbers give us more in- 
formation about large low density regions of the universe, 
and estimates of general characteristics of spatial matter 
distribution derived from such analysis should be essentially 
corrected. 

Results obtained above illustrate the potential and limi- 
tations of this approach rather then give the actual estimates 
of density variations. 



7 SUMMARY AND DISCUSSION. 

In this paper the observed parameters of Ly-a lines are an- 
alyzed and interpreted in the framework of the simple self- 
consistent theoretical model of the DM structure evolution 
(DD99). It is shown that the observed evolution of absorbers 
is sensitive to several random factors and probably does not 
trace directly the evolution of DM component. Our results 
show however that some essential observed characteristics of 
absorbers can be reasonably described even by the statistical 
model considered here. 

The main results of our analysis can be summarized as 
follows: 

(i) The basic observed properties of Ly-a absorbers are 
satisfactorily described by the discussed statistical model of 
DM confined structure elements. 

(ii) The observed characteristics of Doppler parameter, b, 
could be linked to the column density of accompanied DM 
structure elements. It allows us to explain the observed dis- 
tribution of Doppler parameter which is consistent with the 
Gaussian distribution of initial perturbations. The measured 
amplitude of perturbations, Tb, as given by (4.4), is consis- 
tent with that is expected for lower density cosmological 
models. 

(iii) The existence of observed galaxy and simulated DM 
walls as well as the theoretical arguments demonstrate that 
the merging of structure elements is one of the most impor- 
tant factors of structure evolution. The observed character- 
istics of entropy, Fs & S, confirm that such a merging can 
be also considered as probably the main factor of absorbers 
evolution at redshifts z > 2. 

(iv) For the main fraction of absorbers, the weak corre- 
lation of column density, Nhi, with Doppler parameter, b, 
can be reasonably explained by the influence of variations of 
entropy of compressed gas. The smooth observed distribu- 
tion of stronger absorbers, with Nhi > 10 14 cm~ 2 , remains 
in question and must be investigated, first of all, with larger 
simulations. 

(v) The redshift distribution of absorbers does not repeat 
the spatial distribution of galaxies and characterizes mainly 
the matter distribution within larger lower density regions. 
This makes it difficult to reconstruct the spatial DM distri- 
bution from the redshift distribution of absorbers. 

The main statistical characteristics of absorbers distri- 
bution discussed above coincide with published estimates 
(see, e.g., Hu et al., 1995; Cristiani 1995, 1996; Kim et al. 
1997). They are roughly consistent with the expected evo- 
lution of DM structure and can be, in principle, used to 
test and to discriminate different models of structure evo- 
lution. At redshifts z > 3.5 the Doppler parameter and the 
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mean linear number density of stronger absorbers are more 
sensitive to the influence of the initial power spectrum of 
perturbations and to the basic parameters of cosmological 
model. 

The characteristics of DM structure elements discussed 
in Sec. 2.2 and 4.1 were tested against simulated DM dis- 
tribution at z=0 (DD99, DMRT, DDMT). The precision 
reached was ~ 10%. Further application of these methods 
to broader set of simulations at high redshifts will help to 
estimate the unknown numerical factors and to improve the 
proposed here description. The discussed model of structure 
evolution is limited and cannot yet describe, for example, 
the structure disruption. It should be improved at least in 
this respect. 

The available database cannot yet provide a precision 
required for the reliable discrimination of models of struc- 
ture formation. Now the observational data are concentrated 
mainly in the narrow range of redshifts 2.2 < z < 3.5 and 
therefore the quantitative description of absorbers evolution 
is difficult. The wider range of observed redshifts is required 
to improve the description and to obtain more detailed and 
reliable information about the structure evolution and phys- 
ical processes at high redshifts. 
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Appendix A 
COBE normalized power spectrum 

A convenient universal normalization of the power spec- 
trum can be obtained using the anisotropy of Microwave 
Background Radiation measured by COBE. The conve- 
nient parametrisation applied to open universe with mat- 
ter density Q. m < 1 and for the spatially fiat universe with 
fl m + Qa = 1 can be taken from Bunn & White (1996). 

For the Harrison - Zel'dovich initial power spectrum 
with the BBKS transfer function we have: 



a p (z) = 3r (z) — — 



(A.l) 



where a usual definition of moments of power spectrum is 
used: 



x i+ "T 2 (x)dx. 



(A2) 



The function t(z) — toB(z) depends on the model of the 
universe and can be expressed trough Q. m & ft as follows: 



B(z) 



Qm + Oa = 1 



i - n m + 2.2041 + z) 



31 -!/» 



1 + 1.2fi n 



(1 + z)t(z) w 2.73/i 2 fi„ 21 [0.45ft™ 1 + 0.55] 1/3 , z > 1, (A3) 
and for Q, m = 0.3, h — 0.6 we have 

(1 + z)t(z) ^0.29. 
For the open universe with Qa = 0, Q. m < 1 

2.5ft 1 1 



B-\z)^l + 



1 + 1.50, 



-z, 



(1+z)t(z) « l.lh 2 ft^°' 19innm (l + 1.5ft m ), z » 1 (A4) 

and, for example, for ft m = 0.5, h = 0.6 we have 

(1 + z)t(z) w 0.4. 

More details can be found in DD99 and DDMT. 

Appendix B 

Bulk heating and cooling of gas 

The thermal evolution of gas is described by the well 
known equation 

3 dT _ 1 dn b e 7 r 7 n gJ - n 2 b (P fh + f3 ff ) 
2T dt ~ n b dt n b k B T ' [ ' ' 

where e 7 is the mean energy injected at a photoionization 
and radiative cooling coefficients for free-bound and free- 
free emissions can be taken as 

P fb w 3.7 10 _25 (T/10 4 7v-) 1/2 er ff cm 3 /s, 

/3f f w 2 • 10~ 25 (T/10 4 .fO 1/2 erp cm 3 /s. 
If the ionization equilibrium is reached and 

2 ^ 

Urecn b = nmi- 7 
then equation (B.l) can be rewritten more transparently as 



3 dT 



13s = 



2Tn b dt n\ dt 
Pfb + 



J^+ft^/T-l), 



(5.2) 



k B T 



4.2 • 10- 13 {T/W 4 Ky 1/2 cm 3 s- 1 



and T 7 = (5 — 10) • W 4 K for the suitable spectrum of UV 
radiation (Black, 1981). T 7 depends on the spectrum of local 
radiation and can vary from point to point. This heating 
becomes negligible for the cold high density clouds when 
the radiative cooling due to the line emission is essential. 
For the main population of observed absorbers 



bob 



25 - 30km/s, T ob3 w (4 - 6) ■ 10 K < T 7 , 
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and the influence of bulk heating could be essential. 
Eq. (B.2) can be suitably rewritten as follows: 



that allows us to integrate it directly. Using the definition 
of entropy function, F$, as given by (2.20) we obtain: 



where Zf is the redshift of pancake formation. This result 
is obtained under condition of ionization equilibrium. In the 
simplest case of fixed temperatures, T(z) = const < T 7 (z) = 
const, we have 



Fl /2 (z) = Fl /2 (z f ) + a a H [t(z)-t(z f )]\^, (5.5) 



that means the slow growth of entropy and drop of the den- 
sity of absorbers and Nhi- 

The influence of the bulk heating can be enhanced by 
the pancake disruption due to the clustering of both DM 
and baryonic components. This process is usually accompa- 
nied by the adiabatic reduction of temperature in expanded 
regions that accelerates the bulk heating. In this case, how- 
ever, the parameters b and Nhi are strongly correlated. This 
means that the possible contribution of such heating is re- 
stricted at least for richer absorbers. 




T 3/2 /3 s (T 7 /T- 1), 



(5.3) 




= i.5/i _1 c 3 e 6 , 



(BA) 
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